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Abstract. 

The odderon intercept is calculated directly using the wave function recently con- 
structed by R.A. Janik and J.Wosiek. The results confirm their reported value. It is also 
found that their solution for ^3 = does not satisfy the Bose-symmetry requirements. 
Introduction of terms into the trial wave function with an asymptotical behavior simi- 
lar to the Janik- Wosiek wave functions does not seem to improve variational estimates 
significantly. The diffusion parameter is found to be of the order 0.6. 
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1 Introduction. Basic equations. 

1. Recently R.A.Janik and J. Wosiek (JW) published a report on the solution of 
the odderon problem in the perturbative QCD [1]. They derived the odderon intercept 
from the spectrum of the integral of motion q 3 introduced by L.N.Lipatov for the three- 
Reggeon system [2], combined with their earlier solution of the appropriate Baxter 
equation [3]. If one relates the odderon intercept to the odderon "energy" per pair of 
Reggeons e as 

a o (0) = 1 - (3a s /27r)e (1) 
the result of JW for the ground state is 

e = 0.16478... (2) 

Thus they confirm our old conclusion that the odderon intercept lies below unity [3]. 
Their exact value for it is somewhat higher than obtained in variational calculations, 
which gave larger value for e: 0.29 [4] and 0.223 [5,6]. 

In view of a somewhat indirect way of finding the odderon energy, JW suggested 
verifying their result by means of direct calculations of the energy with their wave 
function, using the technical machinery developed in the variational approach. Also, 
to prove that their result indeed gives the minimal value for the energy, it is desirable to 
repeat variational calculation using more general trial functions of the form suggested 
by the JW solution. This short note presents results of these calculations. 

2. The odderon energy can be sought as a ratio 

e = E/D (3) 

where E and D are energy and normalization functionals, quadratic in the odderon 
wave function Z(r,<f>) [5]. Explicitly 

E= £ / dve n (v)\a n (v)\ 2 (4) 

Here 

e» = 2 Re + w) ~ 1^(1)) (5) 

a n {y) is esentially a double Fourier transform of Z(r, 0): 

roo r27T 

a n (v) = / drr- 2 - 2w / d^ 1 ^ 
Jo Jo 
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( n + 1 : . , n — l.\ . n — 1 . . ~ n — 1 . _ . 

+ — + re**(/i - ii/ - — )J (ii/ - + w ~ —) Z ( r ' <t>) ( 6 ) 

where h and /i are the two conformal weights, which are taken to be equal to 1/2 for 
the ground state. D is obtained substituting e n {y) by unity in (1). 

The odderon wave function Z can be considered as a function of z = rexpicp and 
its conjugate. It has a form 

Z(z,z*) = \z(l -^)| 1/3 $(2,^) (7) 
Bose symmetry requires invariance of $ under substitutions 

z-*l-z, z^l/z (8) 

Function $ was explicitly constructed by JW (in the form of a power series). 

Once $ is known, calculation of e reduces to, first, performing the double Fourier 
transform (6) and, second, doing the integration over v and summation over n in (4). 

All the difficulty in the calculation resides at the first stage. Actually one has to 
go to quite high values of n and v to get reasonable accuracy. The Fourier transform 
to such high values of n and v requires much effort. It also requires a very precise 
knowledge of the wave function $ 



2 Odderon ground state. 

1. Because of the latter requirement our first step was to obtain the odderon wave 
function with a higher precision than reported in [1]. To this end we set up a program 
which essentially repeats the procedure employed in [1] and allows to reduce discon- 
tinuities of $ calculated in different variables to values of the order 1CT 9 . The only 
difference with [1] is that we use the basic functions u iy % = 1,2, 3, multiplied by certain 
factors to make them real at points where the transfer matrices are calculated. This 
substantially facilitates achieving the desired accuracy. 

At this step for the integral of motion q 3 and wave function parameters a, f3 and 7 
we obtained the following values (precision 1CT 9 ) 

iq 3 = 0.205257506, a = 0.709605410, (3 = -0.689380668, 7 = 0.145651837 (9) 

However even with these high precision values direct calculation of $ in one of 
the three sets of variables z, 1 - 1/z or 1/(1 - z) fails in the vicinity of the point 
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zq = cxp m I 3 where none of the series converges absolutely. In spite of the fact that zq 
is only a single point in the r, plane, this makes it practically impossible to calculate 
the double Fourier transform for \n\ > 10 and/or \v\ > 5. To overcome this difficulty 
we had to redevelop the function $ around the intermediate point (l/2)z . With this 
redevelopment 100 terms in the series proved to be sufficient to obtain reliable results. 

The Fourier transform itself was performed by interpolating Z quadratically on the 
r, 4> grid and doing the integrals in r and <p analytically. Reasonable results are obtained 
already with a 160 x 160 grid. We however also used 320 x 320 and 640 x 640 grids to 
analyse the precision achieved at this step. 

Even with a 640 x 640 grid the numerical Fourier transform becomes unreliable 
for | n | > 30 and/or \v\ > 15. For such high values of \n\ and \u\ we used asymptotic 
formulas for the Fourier transform, which can easily be obtained from the expansion 
of $ around z — 1. Our final cutoffs were chosen to be |n| < 300 and \u\ < 150., which 
proved to be quite sufficient for the determination of e with a precision 0.001 

2. Results of our numerical calculation of D and E in the region \n\ < 30 and 
\u\ < 15 using an N x iV grid in the r, plane are presented in the Table 1. for 
N = 160, 320 and 640 

Table 1. D and E for the ground state 



N 


D 


E 


160 


1.642162 


0.255480 


320 


1.642085 


0.254852 


640 


1.642056 


0.254632 



To these values one has to add the contributions from the asymptotic region de- 
scribed in the preceding section. They are 

AD = 0.002398, AE = 0.018451 (10) 

Taking the results at N — 640 as the most accurate ones we finally have 

D = 1.644454, E = 0.273083, e = 0.1660 (11) 

Thus our result for e coincides with the value found in [1] up to 0.001. With all 
the difficulties involved in the numerical calculations we consider this agreement quite 
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satisfactory. So direct calculation of the odderon energy confirms the result found by 
JW. 

3. Our previous variational calculations gave very small change in energy as more 
analytic terms were added. To analyse the reason of the about 30% improvement 
given by the JW wave function, we tried to extend our variational calculations to more 
general trial functions as compared to [5,6], whose form is suggested by the latter 
function. At 7^ it possesses an asymptotical behavior at z — > of the same sort as 
the trial wave function introduced by P.Gauron, L.N.Lipatov and B.Nicolescu in [5], 
except of a term which behaves as r 5 / 3 cos 20. To investigate the importance of this 
behaviour we introduced a term into the trial wave function 

ca(z,z*)- 1/6 b(z,z*) 2 (12) 

Here a is the argument in the trial wave function of [5]: 

a = x/y, x= M 2 N 2 , y = (l + \z\ 2 )(l + \z 1 \ 2 )(\z\ 2 + \z 1 \ 2 ) (13) 

where z\ — 1 — z. The argument b is 

b = W /y, w=(l- \z\ 2 )(l - N 2 )(M 2 - N 2 ) (14) 

It is invariant under z — > 1 — z and changes sign under z — > l/z, so that b is invariant. 

At q 3 = a term appears in the solution $ which is unique in the r, <f> plane and 
blows up as r -1 / 3 at r — > 0. Inspection of the functionals D and E shows however that 
it is admissible in spite of the apparent singularity at small r in (6). As mentioned, we 
could not satisfy the Bose symmetry requirements with this solution. Nevertheless we 
tried to estimate its possible significance and so included a term proportional to a -1 / 6 
into the trial function. 

Our final trial function thus included 5 terms, three old ones of the same form as 
in [5,6] and two new ones described above. 

Calculations with this generalized trial function showed first of all that the term 
with A" 1 / 6 presents difficulties for numerical integration, of the same sort that we 
encountered in studying the excited odderon, only much worse. In fact, with this term 
added, we could only calculate the double Fourier transform reliably for |n| < 8, \v\ < 5. 
On the other hand, the term with b lead to no difficulties whatsoever. However in both 
cases one finds that the new added terms bear no influence on the value of energy. At 
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the minimum value for the functional E the coefficients before them turn out to be 
quite small and the value itself is only a few percent lower than without the new terms. 
So our conclusion is that simple addition of new terms into the trial wave functions even 
with two independent arguments a and b described above does not improve variational 
estimates. It is a factorized form of the Janik-Wosiek wave function which allows to 
make energy substantially lower. These calculations also make us believe that the JW 
function indeed belongs to the odderon ground state. 

3 Other eigenvalues of q 3 

1. We have also tried to check the result of [1] for the excited state with the next higher 
value of iq%. Unfortunately in this case calculations proved to be still more difficult 
and we could not arrive at a result of a convincing accuracy. 

Our precise calculations of the wave function gave for this state: 

iq 3 = 2.343921063, a = 0.391855163, (3 = -0.0533712012, 7 = 0.918477570 (15) 

Numerical calculation of D and E in the region |n| < 15, \u\ < 15. gave results 
presented in Table 2. 

Table 2. D and E for the state with iq 3 = 2.34... 



N 


D 


E 


160 


2.92863 


6.08989 


320 


2.80693 


5.27381 


640 


3.05939 


6.96764 



As one observes the achieved accuracy does not exceed 15%. Analysing these num- 
bers one can see that all the error comes from the region of maximal \n\ and \v\ where 
the double Fourier transform is apparently performed inaccurately. /,From these num- 
bers we can only conclude that for this excited state 

e~2. ±0.3 (16) 

In [1] the found value is 1.71231... Our result does not contradict this number. 

2. We have also studied a possible solution for q 3 = 0, reported in [1]. However in 
this case we were not able to construct a wave function unique in the r, <fi plane and 
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satisfying the necessary symmetry requirements. If, following [1] we seek $ in the form 

$(z,z*)=uAu (17) 

where Ui(z) are the three basic solutions of the eigenvalue equation for q 3 and require the 
matrix A to have certain matrix elements equal to zero to make $ finite, then at ^3 = 
we find it impossible for A to be invariant under the basic symmetry transformations 
of z. This seems to be related to the fact that at ^3 = one may also take A 33 7^ 0. In 
any case there does not seem to exist solutions for A which preserve the required Bose 
symmetry. 

Thus our conclusion is that the g 3 = state reported in [1] does not correspond to 
any physical odderon state. 



4 "Moving" odderon 

For conformal weights h = \ + ia the odderon energy is supposed to behave at small 
a as 

e(a) = e + aa 2 (18) 

where eo is the value (2) and a is a parameter which determines diffusion of the odderon 
wave function in the momentum space. This parameter has been long known for the 
pomeron to be 14£(3) (in units 3a s /7r). It is of certain interest to find a for the odderon. 

To this aim we first found the parameters of the odderon wave function for various 
(small) a using the same method as employed for h — 1/2. Our results are presented 
in Table 3. The value of iq 3 turned out to be real for arbitrary a, whereas, with a 
chosen to be real, both (3 and 7 result complex. We chose a — 1 



Table 3. Odderon parameters for h = \ + ia 



a 


iq3 


P 


7 


0.01 


0.205306079 


-0.971740164-i0.014404102 


0.205305637-i0.000425478 


0.1 


0.210089247 


-0.995153863-10.142974530 


0.210052319-10.003938872 


0.2 


0.224303315 


-1.060790013-i0.281079150 


0.224222881-i0.006006327 


0.3 


0.247227544 


-1.156524786-10.415163678 


0.247186043-i0.004529717 


0.5 


0.316528176 


-1.395571390-i0.695891904 


0.316214188+i0.014095104 


1.0 


0.619239545 


-2.044631201-il.672240784 


0.591391973+i0.183611401 
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Inspecting these figures one immediately notes that |7| = iq%. This relation was 
predicted (for real 7) by L.N.Lipatov [7]. 

With the odderon parameters found we calculated the odderon energies directly, 
using the same techique as for a = 0. With a different from zero calculation becomes 
still more cumbersome and time and memory consuming due to lack of certain sym- 
metries and overall complex arithmetics. For these reasons we had to limit ourselves 
with a maximal 160 x 160 grid in the r, plane and neglected the contribution from 
the asymptotical region n > 30 \v\ > 15. Our results are shown in Table 4 together 
with the ones just obtained via the solution of the Baxter equation [8] 

Table 4. Odderon energies for h — \ + ia 



a 


e 


6 [8] 


0.0 


0.1534 


0.16478 


0.1 


0.1597 




0.3 


0.2085 


0.21777 


0.5 


0.2980 


0.30523 


1.0 


0.6269 


0.63228 



Our energies lie a little below the ones obtained from the Baxter equation, which is 
natural since we have neglected the asymptotic part of the n, v region in (4). Having 
this in mind we find a complete agreement between our direct calculation results and 
the ones based on the Baxter equation. 

/From our energies we find for the parameter a in (18) 

a = 0.61 

More precise energies found in [8] lead to 

a = 0.605 

Note however that already at a — 1 the approximation (18) breaks down and more 
powers of a 2 are needed to describe the energy behaviour. It is interesting that the 
parameter a for the odderon is much smaller than the one for the pomeron. In fact 
their ratio is of the same order as the ratio of corresponding energies. 
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